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In extra-dimensional scenarios oscillations between active and sterile neutrinos can be governed 
by a new resonance in the oscillation amplitude. This resonance results when cancelation occurs 
between two phase differences, the usual kinematic one coming from the neutrino mass-squared 
difference, and a geometric one coming from the difference in travel times of the sterile neutrino 
through the bulk relative to the active neutrino confined to the brane. In this work we introduce a 
specific metric for the brane-bulk system, from which we explicitly derive extra-dimensional geodesies 
for the sterile neutrino, and ultimately the oscillation probability of the active-sterile two-state 
system. We find that for an asymmetrically-warped metric, the resonance condition involves both the 
neutrino energy E and the travel distance L on the brane. In other words, the resonant energy may 
be viewed as baseline-dependent. We show that to a good approximation, the resonance condition 
is not on E or on L, but rather on the product LE. The model is rich in implications, including the 
possibility of multiple solutions to the resonance condition, with ramifications for existing data sets, 
e.g., LSND and MiniBooNE. Some phenomenology with these brane-bulk resonances is discussed. 

PACS numbers: 13.15.-|-g, 14.60.Pq, 14.60.St 



I. INTRODUCTION 

In theories with large extra dimensions, the Standard Model particles are typically confined to the 3-1- 1-dimensional 
brane, which is embedded in an extra-dimensional bulk [1, 2]. Singlets under the gauge group however are allowed to 
travel freely in the bulk as well as on the brane. 

It has been argued that the LSND neutrino oscillation anomaly [3] and the MiniBooNE [4] null result [23] might 
be explained by a brane-bulk resonance in active-sterile neutrino oscillations [5]. The resonance arises due to the 
additional phase difference 5{Ht) — tSH + H5t induced when the sterile neutrinos take temporal shortcuts through 
an extra dimension. In general, _ff is an n x n matrix for n neutrino flavors. 

Thus, in these models there are two sources of phase difference, the standard one t5H = L/S.m? /2E, and a new 
one Ht{St/t) arising from temporal shortcuts through the bulk available to gauge singlet quanta. The two phase 
differences may beat against one another to produce resonant phenomena. Here, as in [5], we will focus on the gauge 
singlet sterile neutrino [6]. 

For simplicity, the shortcut parameter e = 6t/t in [5] was taken to be a constant. It was given a physical interpre- 
tation as the "aspect ratio" of a typical brane fluctuation; quantitatively, it was related to the fluctuation geometry 
as e = (^1^) , where Su is the amplitude of the brane fluctuation in the bulk direction u, and Sx is the length of the 
fluctuation along the brane direction x. It is the purpose of this article to extend the idea in [5] to an asymmetrically- 
warped metric framework for the extra dimension. Significant new physics emerges from this generalization. It will 
be seen that the shortcut parameter is no longer a constant, but rather depends on the neutrino propagation distance, 
i.e., on the experimental baseline length L. Larger baselines allow more effective shortcuts (larger e), which leads 
to correspondingly smaller resonance energies. This effect can be traced back to the travel times of the neutrinos. 
Longer neutrino travel times associated with longer baselines on the brane allow the off-brane geodesic of the sterile 
neutrino to plunge deeper into the bulk and experience a greater warp factor. The smaller resonance energies, in 
turn, suppress active-sterile neutrino mixing at large energies. Thus, one expects little or no active-sterile oscillations 
at the relatively large baselines accompanying atmospheric and astrophysical neutrino sources. 
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For the asymmctrically-warpcd metric which we investigate, a remarkable feature emerges. It turns out that to a 
good approximation, the location of the resonance depends on the product of baseline and energy, LE = (LE)-^^^, 
where {LE)-[ics is a constant involving the warp factor and neutrino mass and mixing angle. From this result one 
infers that experiments with LE well below (Li?)Res should observe a 5D version of active-sterile vacuum oscillations 
(presented in section VB), while experiments with LE well above (Li?) Res should have active-sterile mixing suppressed 
below observable levels. 

Another feature of the new phenomenology presented herein, is the possibility of multiple brane-bulk resonances, 
with interference among them. This is because there are multiple extra-dimensional paths contributing to the prop- 
agation amplitude of the sterile neutrino, each with a different intrinsic resonant energy. The result of a quantum 
mechanical sum of paths is a richer picture for the resonant behavior than that of the sharp resonant energy in [5] , 
or than that which occurs in neutrino propagation in a nonzero matter density [8, 9]. 

In the next section we review the physics of the brane-bulk resonance, to set the stage for the further developments 
presented in this work. 



II. TWO-STATE ACTIVE-STERILE OSCILLATIONS IN A BRANE-BULK SCENARIO 

Following [5], we consider two mass eigenstates i^i and and a small active-sterile mixing angle 0. The active-sterile 
mixing angle 9 relates the flavor eigenstates Va, fs with the mass eigenstates v\, V2 via a unitary transformation 

2 

I fa) = ^Z7*j| fi) with a = a, s (1) 

and 

U=( ^^"^.V (2) 

\^ — sm f cos J ^ ' 

With small mixing angle, the state Va is mostly Vx, and the state Vg is mostly y2- 

The Schrodinger equation in flavor space describing two-state neutrino oscillations reads 

dt V '^s{t) [ >ys{t) J ' ^"^^ 

where the relevant Hamiltonian is given by 

Am^ f-cos2e sin 26* 

The "induced mass-squared difference" is 

B{E) = 2E^e, (5) 

and we have dropped the irrelevant piece of H equal to E + {m^ -\- mD/AE times the identity matrix. The parameter 
e = dt/t is the shortcut parameter, which parameterizes the relative difference in travel times of active and sterile 
neutrino flavors. It will be derived explicitly for an asymmetrically-warped metric in section III. The novel feature 
of this work will be that e, and the resonant energy E^Res that results from it, will show dependencies on the baseline 
length L. 

This Hamiltonian system bears much resemblance to the Mikheev-Smirnov-Wolfenstein (MSW) Hamiltonian [8, 9] 
describing neutrino oscillations in matter. Accordingly, we shall henceforth refer to 6 and Am^ as the 4D vacuum 
mixing angle and 4D vacuum mass-squared difference, respectively. In the brane-bulk system, the shortcut in the 
extra dimension has been parameterized by an effective potential with nonzero sterile-sterile component in the ffavor 
space Hamiltonian. This term is the analog of the Vg-t^e potential term in the MSW Hamiltonian, induced by coherent 
elastic forward scattering of neutrinos on electrons. 

Note, however, that there are three important differences between the brane-bulk system and the matter system. 
First of all, as previously mentioned, gravitationally determined geodesies do not discriminate between particle and 
antiparticle and the effective potential is therefore the same for neutrino and antincutrino. Secondly, the energy 
dependence is stronger for the brane-bulk system than for the matter system. In the brane-bulk system, the induced 
mass-squared difference varies as E^, whereas in the MSW system [9] the variation is proportional to E. The third 
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difference is that in the brane-bulk system, there is no time or space dependence in the effective potential, and 
therefore none in the Hamiltonian. 

Diagonahzation of the 2x2 effective Hamiltonian proceeds analogously to the MSW system. One finds a new 
mixing angle 9, expressible in terms of 4D vacuum values as 



tan 261 = 



tan26» 



1 - 



(6) 



with 



Am^ cos 26 



As with MSW, neutrino mixing can be maximal {9 
Eq. (7) can be rewritten as 



2e 



7r/4) even for a small 4D vacuum mixing angle. 



(7) 



e = 



cos 29 Am? 



(8) 



which shows that e will be very small for E^^^ ^ Arr? . In [5], Arr? was taken to be ~ eV^ to accommodate the 
LSND data, and -Ercs was taken above 30 MeV and up to the MiniBooNE excess region of ~ 300 MeV. Consequently, 
e turned out to be about 0.5 x lO^^^"*"^. 

The existence of a resonant energy naturally divides the energy domain into three regions. Above the resonance at 
E ^ £^Res, the sterile state decouples from the active state as sin^ 29 0, and therefore oscillations are suppressed 
in this regime. At resonance, i.e for E = E'ros, the mixing angle 9 attains a maximum. Finally, below the resonant 
energy, the oscillation parameters reduce to their 4D vacuum values. 

The new eigenvalue difference 5H is given by 



A ^2 / 

dH = —— \ / sin^ 261 + cos2 26* 
2E v 



-^Ros 



and one obtains the usual expression for the flavor-oscillation probability 



. 2^;^ . 2L6H 
Pas = sin^ 26* sin^ — — 



(9) 



(10) 



with the expression 



sin^ 29 



sin^ 29 



sin^ 26* + cos2 26* 



1 



(11) 



following from Eq. (6). 

Such is the picture for a sterile neutrino traversing a unique shortcut through the bulk. In the next section, we avoid 
the nebulous feature of a "typical fluctuation" of the brane in the bulk by introducing a metric for the brane-bulk 
system. This will allow us to solve explicitly for the brane-bulk geodesies (plural) of the sterile neutrino. We will find 
that there are multiple geodesic paths available to the sterile neutrino. Each path corresponds to an amplitude, and 
the multiple paths necessitate a quantum mechanical sum, which complicates but enriches the analysis. 



III. A METRIC AND GEODESICS FOR STERILE NEUTRINO SHORTCUTS 



A sterile neutrino which propagates in the bulk as well as on the brane will in general have a different trajectory 

compared to the active flavor which is confined to the brane. If the extra dimension is "warped" relative to the brane 
dimensions, then the geodesic of the sterile neutrino will be shorter than the geodesic of the active neutrino. For an 
observer on the brane, it will appear as if the sterile neutrino has taken a shortcut through the extra dimension. Such 
apparent superluminal behavior for gauge-singlet quanta has been discussed earlier for the graviton [10, 11]. 
Specifically, we consider an asymmetrically warped 4 -I- 1-dimensional spacetime metric of the form 

3 

dr^ = de-Y^ V^{u) {dx')^ - du", (12) 

i=l 
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where u is the extra dimension and 'q{u) is the warp factor. Since space but not time is warped, the warp is termed 
"asymmetric" [11, 12]. The Standard Model neutrinos Uve on the brane while the sterile neutrinos propagate freely in 
the extra dimension. Our observable brane is the u = submanifold, which wc take to be globally Minkowskian, i.e. 
ri{u = 0) = 1. We choose a warp factor of the form r]{u) = e"*^'"' , with k assumed to be a (presently) unknown constant 
of Nature with dimension of inverse length [24]. Standard Model neutrinos live in the 4D Minkowski spacetime, while 
the sterile neutrinos experience the full five dimensional metric. We may choose the direction of the brane component 
of neutrino velocity to be along x; this allows us to set dy and dz to zero from here on. Our line element is reduced to 

dT^ = de - e-2'=l«l dx^ - du^ . (13) 

We choose the absolute value of u in the warp factor such that the bulk acts like a converging lens on both sides of 
the brane, returning the sterile neutrino to the brane where it may interfere with the active neutrino. 

Our task is to calculate the difference 5t = te'm the propagation time between the Standard Model (active) neutrinos 
which propagate only along our brane, and the sterile neutrinos which take shortcut excursions through the warped 
extra dimension. For the active neutrinos, we have immediately the Minkowski space result £°^^'^'^ = Ljji, where /? is 
the speed of the neutrino. For the sterile neutrinos, we must find and solve geodesic equations. It is well-known that 
geodesic equations may be found by treating dr^, or equivalently, [dr j d.Xf' , with A any afRne parameter which marks 
progress along the geodesic, as a Lagrangian subject to the extremization conditions given by the Euler-Lagrange 
equation. Here we have 

There result three Euler-Lagrange/geodesic equations, one for each of the variables {t^x,u\. They are the following 
set of coupled nonlinear differential equations 



= (15) 



Afe--H^) =0 (16) 



d\ \ dX 

The ± signs refer to trajectories in the positive or negative half-planes of u. By integrating Eq. (15) twice, one readily 
infers that A and t are proportional to one another. We call the proportionality constant dt/dX = 7. It is most easily 
obtained from the line element in Eq. (13), by taking dr to be dX, as is appropriate for massive particles. One gets 

and we have denoted some occurrences of d/dt by an overdot, to streamline notation. Initially, the neutrinos are 
created on the brane, where m = 0, so we have 



7^ = ^^, with (3=^x1+ ill. (19) 

A variable with a zero subscript denotes an initial {t = 0) value. The interpretation of 7 and (3 is now clear: /3 is the 
neutrino speed, and 7 is the usual boost factor, equal to E/m. The only difference for the sterile neutrino versus the 
active neutrino is that the sterile is allowed to have an initial velocity component Uq in the bulk direction, while the 
active neutrino cannot. 

Translational invariance is maintained on the brane and in fact, the Minkowski metric on the brane assures that 
Lorentz invariance is maintained on the brane. Therefore, momentum components are conserved on the brane, and 
we cannot generate a nonzero ii on the brane except as an initial condition. Thus, we must have a nonzero iio in order 
for the sterile neutrino to leave the brane and enjoy the bulk. The uncertainty principle applied to the u dimension 
allows for such a nonzero velocity. In fact, it necessitates some uncertainty in uq. With a brane thickness of Aw,, one 
expects an initial Pu = m'yuo ~ 1/ Au. The path integral approach to quantum mechanics, a topic to which we will 
return later in this paper, gives a similar result. In the path integral approach, all uq are allowed, with each value 
weighted by the free-particle action {S = m J dT)/h. The dominant paths are the classical cxtrcmac of the action, 
and their width in action, of order h, incorporates the quantum mechanical uncertainty. Momentum conservation in 
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the ?i-direction is a non-issue, as translational invariance in the w-direction is broken by the brane itself; alternatively, 
one may simply view the recoil of the brane as compensating any momentum change in p„. 

Since the set of geodesic equations are homogeneous in A, we may therefore replace A with t, and we do so. Thus 
we arrive at two coupled geodesic equations, 



dt 



-2k\u\ 



dt J 



dt^ 



dt 



= 



0. 



(20) 
(21) 



For Eq. (20) to be satisfied, the expression in parenthesis must be a constant, independent of the spacetime variables. 
By taking u to zero, it is clear that the constant is xq, the initial velocity component of the sterile neutrino along the 
brane. So we have 



at 

Using this result, the geodesic Eq. (21) can be rewritten as 

f^±kxle''^\-\=0. 



(22) 



(23) 



This has a standard form for a second order nonlinear differential equation, which can be once-integrated and recast 

as 



T^kxl 



u{t) 



du e^'^l"! . 



(24) 



The trajectory of the sterile neutrinos is contained in Eqs. (22) and (24). From these equations, we may form the 
ratio x/ii = dx/du and solve for x{u) or u{x). We obtain 



u{x) = ±— In fl + k^x{l - x)] 
2k 



(25) 



where / = 2\uo\/kxo is the distance along the brane at which the sterile neutrino trajectory returns to the brane. For 
small kl <^ 1, this geodesic solution is the particularly simple parabola, u{x) = ±| x{l — x). 

Notice that the dimensionful warp parameter k sets the scale for all length measurements. Accordingly, Eq. (25) is 
especially simple when lengths are measured in units of k~^. 

The geodesic is symmetric about the midpoint at x = 1/2, which implies that the angle of intersection of the sterile 
neutrino trajectory coming back to the brane is the same as the initial exit angle, given by arctan(uo/io)- From this 
we infer two features of the geodesic. The first is that the exit angle is also equal to arctan(A;?/2). This presents 
a physical picture for the strength of the warp parameter k. The second feature follows from the symmetry of the 
geodesic and the "left-right ?i-symmetry" of the warp factor. It is that \u{x)\ as given in Eq. (25) repeats itself 
indefinitely, alternately on opposite sides of the brane. Thus, the geodesic is periodic, with I specifying the length of 
a half-cycle. A half-cycle of the geodesic ku{kx) is pictured in Fig. 1 for some choices of kl. The depth of penetration 
into the bulk is 



Iwrnaxl = \u{x = l/2)\ = ^ In 



1+ ^ 



(26) 



This result provides another picture for the strength of the warp parameter. At small kl we have fcumax ~ {kl)'^/8, 
and at large kl we have fcwmax ~ ln(fc?/2). 

It will be useful to have expressions for the speeds iio and xq in terms of the geometric quantities k and I. From 
I = 2uo/kxo above, and the constant of motion = 11^ + Xq, one easily finds 



Xo 



f3 



yr7(ff ' 



and 



Uo 



/3(f) 



(27) 



It is important to keep in mind that the distance I as measured on the brane is a function of the initial values of the 

velocity components along the extra dimension and along the brane. Further in the paper we will mostly refer to 
dependencies on I with the implicit dependence of I on the initial conditions. 
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The travel time t^^^^ of the sterile state in the bulk is determined by substituting the expressions for x and u{x) 
given in Eqs. (22) and (25) into 



dx 



dx 



x{u{x)) Jq xo[1 + k^x {I — x) 



kl f l \ f 2 kl 
arcsmh — = \ — \ — arcsmh — 
,2 2 \kl 2 



(28) 



iofcVl + (f)^ 

The travel time for the Standard Model neutrinos confined to the brane is 

^branc(-;) = ///?. Thus, the Sterile neutrino 
appears to a brane-bound observer to travel super luminally, exceeding light-speed by the factor (^ arcsinh^)- One 
can now write down the "shortcut parameter" e{l) in terms of the brane distance I as 



-ibranc xbulk / R\ 

eil) ^ i = 1 - ( ^ ) t*^""^ = 1 

^ ' y- brane \ / 



^) arcsinh 

kl V 2 



(29) 



An important feature has emerged. The shortcut parameter depends on the length / of the geodesic as seen from the 
brane. According to Eq. (26), the geodesic for large I dives deeper into the bulk, where the warp factor exponentially 
decreases the travel time. The dependence of e on kl is shown in Fig. 2. Focus here on the curve labeled n = 1\ the 
other curves will be explained shortly. 

Phenomenologically, we are interested in whether the sterile neutrino returns to the brane at a fixed baseline L, 
where an experiment might be located. If the sterile neutrino returns to the brane at L, then its wave packet will 
interfere with the active neutrino packet to produce flavor oscillations. If the sterile neutrino does not return to the 
brane at L, then there is no flavor interference. The baseline L and the return length I are related in a simple way. 
Given the periodicity of the geodesic, a sterile neutrino will interfere on the brane at L if L is a multiple of the 
half-cycle length I, i.e., if 



n 



(30) 



where n is any positive integer. The travel time and therefore the shortcut parameter for a sterile neutrino which 
enters the detector upon intersecting the brane an v}^ time is simply obtained by replacing I with L/n in Eqs. (28) 
and (29), respectively. The results are 



t 



bulk 



{L) 



Here we have defined the dimensionless "scaling variable" 



and 



e„(i) 



( — J arcsinh ( — 
(- 



L 



(1 - en) , 



arcsinh 



kL 



(31) 
(32) 

(33) 



This variable will appear throughout subsequent equations. Its dimensionless feature proves useful when a perturbation 
series in v or v/n is desired. Fig. 2 shows the dependence of e„(i) on v = kL/2 and n. It is seen that e grows 



7 




V = kL/2 



FIG. 2: The shortcut parameter as a function of scaled baseline length v = kL/2. Curves are parametrized by geodesic mode 
number n = 1, 2, 5, 10 (from top to bottom). 



monotonically with v, and for fixed w, shrinks monotonically with geodesic mode number n. Since we have seen in 
Eq. (7) that the resonance energy -Ercs varies as l/^/e, the energy i?Rcs will shrink monotonically with increasing L, 
and grow monotonically with increasing n. 

Anticipating the next sections of this paper, we expand Eq. (32) in v/n. The result is 



1 



9 /vy 
20 \n) 



(34) 



The leading behavior tells us that small e and small are synonymous, and that e„ « (kL/n)^ oc L^, for small 

v/n. Therefore, for small v/n one has E^cs{n) ~ 2nV3Am'^ cos29/kL as the resonant energy of the n*'' mode. 
Also useful will be the generalization of Eq. (27), obtained by setting I — L/n. The generalization is 



and 



HI) 



Pv 



(35) 



Eq. (35) reveals that the initial angle between the sterile trajectory and the brane, ao, is given by tanao = juol/io 
v/n. 



IV. SUMMING THE GEODESIC PATHS 



Each geodesic, labeled by integer n's, presents a possible path for the sterile neutrino. According to quantum 
mechanics, we must sum over these possibilities (amplitudes A(n)), and square the total amplitude {A = A(n)). 
In this way, our calculation becomes a semi-classical approach to path-integral quantum mechanics. Each geodesic is 
a classical path, and so is also an extremum of the path integral. 

What is absent at this point is the weight to be associated with each classical path. There are actually three weights 
to be considered. 



A. Path-integral weight 



The first is the usual quantum mechanical weight e " , where 



Sn — m dr 



m 



dt 



m 
~E 

bulk 



.bulk 



m L 

(3E 



(1 - e„) 



(36) 



is the n mode action for a free-particle. The variables i„ and e„ are given in Eqs. (31) and (32) above. The small 
V expansion of Sn is 



Sn = 



(3E 



1 - 



D \n/ 



(37) 
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To leading order, Sn is independent of n. Therefore, for small v/n, the uniniodular weight e''^" ean be omitted 
from the sum. This means that for small v, all amplitudes A(n) are relatively real and so interfere constructively. 
On the other hand, for intermediate and large values of v/n, the weight is n-dependent and must be retained. For 
intermediate and large v/n, interferences among modes depend on the phase 5„. However, from Eq. (58) we infer that 
iJ^gg(n) = Ato^ cos 20/2 e„, so we may infer that for a resonant energy in the interesting range E^gs{n) ^ VAm?, 
we are in the small e„, equivalently, small v/n domain, and the weight e'^" can be omitted. 



B. Distribution of initial uo 



The second weight concerns the a priori amplitude for nonzero values of uq. According to Eq. (35), there is a 
one-to-one correspondence between v/n = kL/2n and uq. Thus, each mode intersecting the brane results from a 
unique initial value of uq. Here we make our most model-dependent assumption, namely that the initial amplitudes in 
momentum pu = mjiio transverse to the brane assume a normal distribution about p„ = 0. The arguments favoring 
this choice are standard: Ground state wave functions often have a GauBian distribution about the mean; random 
fluctuations are distributed in Gaufiian form; while initial momenta may prefer the brane direction, the uncertainty 
principle requires a nonzero p„ component as well. Thus, we assume for the uo-distribution the following normalized 
Gaui3ian 

1 -14) 

dNaiPu) = -j=- e i dpu , (38) 



which satisfies the usual normalization condition 

*>+oo 



/ 



1 



2Tra 



dpu = l. (39) 



The new parameter a is the rms value of p„. Heisenberg's uncertainty principle relates a to the uncertainty in the 
initial u-position Au via 

^ > TT^- (40) 

2Au 

We will take a ~ 1/Au, and interpret the uncertainty Am as some measure of the thickness of our brane. Another 
estimate would be to set a equal to the only extra-dimensional scale in the problem, the warp parameter k. If the 
extra dimensions are compactified at a length scale ■Ad'^^ , then certainly Au < ■Ad'^^ , or equivalently, a > Add- 
Reference [13] provides a lower bound on the scale of the extra dimension Md — 10~^ eV, or M'^^ ~ 0.2 mm. While 
(J ~ 1/ Au is a parameter which in principle can be fit to data, e.g. to the height of the resonance oscillation peak, as 
we show later. We will use a fiducial scale of one micron for Au in our numerical work. This translates into a value 
of 0.2 eV for the rms momentum a, which is coincident ally (or not?) of order of the scale of neutrino masses. 
Employing m.'y = E and Eq. (35), the momentum component along u can be written as 

Pu = jmuo = , - „ . (41) 
V?i + v-' 

Substituting for p„ with Eq. (41), our diff'erential Mo-distribution in Eq. (38) can be rewritten as 



y/2na " dn 



z;^ - ■ - / 9 , i-dn). (42) 



1 A^^m^] E(5^ 
Integration of {—dn) is to be interpreted as 

M poo oo 

2 / {-dn) = 2 rfn = 2 V An , (43) 
Joo Jo 
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where the factor of 2 accounts for the negative integers which label the geodesies whose first half-cycle is in the 
negative u half-plane, and An is a measure to be associated with each geodesic mode. Thus, we arrive at 



JdNoin) = 



main) = r-^Y^^r 



[nP- + ^2)3/2 



(44) 



If An were \mity for all n, then the sum amounts to approximating the Gaufiian integral by a histogram, with each 
integer specifying a bin of unit width. This would yield Ylm ^ ~ which means that every neutrino which leaves 
the brane appears at a; = L. Since in general some sterile neutrinos leaving the brane will not arrive aX x = L, this 
choice of measure for An cannot be so simple. The measure for An. our third weight, will be discussed shortly. 

Guesses for the initial distribution of uqS other than the GauBian form presented above are possible. A simple 
example would be to distribute uo uniformly over the interval [— Ug"*^, w™**^]. With this guess, one has the flat 
distribution 



dNpin) 



dNf 



diif] 



dn 



1 



(3vn 



{-dn) . 



diia J \dn J '~ \2u^''''' J \{n^+v^)^/^ 
Employing the left-right bulk symmetry of our metric and moving to discrete notation, this may be written as 



/ 



dNpin) 



.^Ar 

n=l 



(„2 ^ ^2-)3/2 



(45) 



(46) 



Similarities and differences between Eqs. (44) and (46) arc evident. The advantages of the ?io-distribution in Eq. 
(44) are that the Gaufiian form is motivated by physical arguments given earlier, and that the distribution is sensitive 
to the brane thickness Au ~ l/cr. One expects the brane thickness to enter into the physics, and it does so here. 
A disadvantage for Eq. (44) is that there is one more parameter assignment to be made, the value for cr. However, 
since a entered originally as part of the normalization of the continuous Mo-distribution, and therefore vanished upon 
integration, we do not expect much residual sensitivity to a when we move to the discrete weighted sum. The 
advantage of the lio-distribution in Eq. (46) is its sheer simplicity. This distribution may be motivated by appealing 
to path-integral quantum mechanics, where a priori all initial velocities are allowed with equal weight. In fact, we 
may write both distributions dNo and dNp in a similar way: 



n^ + v"^ 



-( V \ 



D = Gaufiian 



D = Flat . 



(47) 



Written in this form, the similarities and differences are more apparent. Both bracketed expressions are bounded 
from above by 0(1). In fact, the slope of the Gaufiian is bounded by 1/Ve ~ 0.7, and both are small for small v. 
The differences appear at very small v and at large v. At very small w, the bracketed expression is simply v/n for 
dNp, while for dNa it is ^J^Jn 0Ev / an, and so is dependent on the brane width Au, ^ l/cr for dNc. At resonance, 
the bracketed expression for dNc is yj&l^m? cos26'/7r (/3/cr) ~ m/a. Also, dNc is small for large v (the tail of the 
Gaufiian), whereas dNp is 0(1) at large v. As mentioned above, for a resonance to occur in an interesting energy 
range requires that v be small, so small v becomes our focus. 

We will analyze the implications of both ug-distributions. For the Flat Uo-distribution, we will set ug"'^ — (5 \n our 
numerical work. 



C. The weight of An 

We come to the third and final weight. An, characterizing the thickness of the classical path in n-space. According 

to path- integral quantum mechanics, the thickness of the path is determined by the variation of e*^/'' about the 
classical extremum e^Sd/h^ ^qy \S — Sci\ > h, variations are rapid and the integration averages to zero. Thus, A^* ~ h 
characterizes the thickness of the path. So we have 
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where S{n) is viewed as a continuous function of n. Non-integer values of n describe geodesies which interpolate 
between the integer-labeled modes, i.e. they describe classical paths which do not intersect x at L. But if the 
dependence of 5 on n is mild, i.e., if dS/dn is small, then the quantum mechanical uncertainty allows the interpolated 
paths to intersect x a.t L, too. Sd is the extremum of the multi-functional S{t;x,u) that characterizes a classical 
path. But Sci is not an extremum of S in the variable n. We assume the dependence of 5 on n is sufBciently mild 
that dS/dn ~ dSd/dn. Prom Eqs. (36) and (32) we easily obtain 



dSci jn) 
dn 



2m^ ( V 
-— "T arcsmh— 
Epk \ n 



(49) 



Of course. An cannot exceed unity, the distance between successive integers. So we arrive at our definition of the 
measure 



An = min < 1 



Lm? 



arcsinh- 



(50) 



The small v expansion of Eq. (50) is 



An = min 



9 /vy 



for t;^ < 1 . 



(51) 



It is useful at this point to introduce the usual vacuum oscillation length for neutrinos on the 4D brane A4D = 
4:-kE/Aw? = 2.48km(i;/GeV)(Am2/eV^)-^ Then the prefactor E/3/Lm^ may be written as /3/47riv{'*^\ assuming 
and where we have introduced the dimensionless 



(4D) 



L 
A4D 



LAm^ 
A-kE ' 



(52) 



which is equal to the number of vacuum oscillations in the baseline interval to L. N^^^'' is a useful variable in that it 

is effectively bounded by 100, at most, if oscillations arc to occur. This bound arises from "experimental dccohcrcnce" 
as follows [25]: The 4D vacuum oscillation phase is ^ = L Am^ /2 E. Therefore, experimental uncertainties in L and 
E introduce an uncertainty in the phase of 




(53) 



But (f) = 2ttN[^"\ so even if the relative experimental uncertainties in L and E were as small as a few percent, the 
phase uncertainty becomes total, i.e. dcf) = 2tt, before N^^^ reaches 100. So by writing the explicit Lm^/E de- 
pendence in terms of the dimensionless n'^^\ we are assured of a pure number whose value cannot exceed 100 if 
coherence is to be maintained. Thus, the option (^u^) ^^^^ greater than ~ 3n^/(4007rw^) if coherence 

is to be maintained. This exceeds unity, the alternate option, if < n'^/400. So for v < if coherence is to 

be maintained, then we have An(n) = 1. If v < 1/20, then An(n) = 1 for all n. The physics of this result is that 
for V < n^/^/20, the action is sufBciently flat in n (Recall that Sn as given in Eq. (37) is nearly independent of n for 
small v/n.) that all neutrinos leaving the brane on trajectories labeled by the interval [n, n -|- 1] return to the brane 
at X = L [26]. 

Eqs. (26) and (33) give the maximum classical penetration depth into the bulk as 



which has a small v/n expansion 



Lv 



2 \n) 



(54) 



(55) 



Thus we see again that for small v, the geodesies never wander farther than vL/A from the brane. In this sense it is 
not surprising that all neutrinos leaving the brane at small angles arctan u u have an amplitude to be found on the 
brane at a; = L. 
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If N^''^^ does exceed 0(100), coherence is lost, oscillations are unobservable, but oscillation-averaged mixing persists 
as the vestigial signature of oscillations. In fact, it is N^^'' which is relevant for determining coherence or decoherence 
in the present model. Setting Asd = 27r/(5i/„ leads to the relation 



N- 



{5D) 



N 



{4D) 



/sin^ 2e + cos2 29 



(56) 



Except near a resonance, A^^^^^-* is a good guide for N^^\ At a resonance, the value of N^^^ overestimates the 
number of oscillations by the factor of 1/ sin 20. We return to the possibility of decoherent mixing in the next 
subsection. 

D. Putting all the pieces together 

The total amplitude for the active-sterile neutrino oscillations is the sum of the amplitudes A{n) for each mode n. 
Each A{n) includes the three weights 6*^"=''"^ dNo, and An. We have 



vn 



n=l 



n=l 



(n2 +^2)3/2 



2 /3E 



sin2^„ sin 



(57) 



with sin20„ and SHn obtained from sin 20 and SH in Eqs. (11) and (9) by replacing .ERes with the resonant energy 
of the mode 



I Am2 cos 26 



Am? cos 26 



2e„ 



2[l-(^)arcsinh(^)] " 



(58) 



The bracketed expression in Eq. (57) pertains to the dNc weighting of the Mo-distribution. For the dNp weighting, 
this bracketed expression is replaced by The probability of oscillation is the square of the amplitude of 

oscillation, 



p _ M |2 _ 
^ as — l-'^asl — 



(59) 



Recall that if TV^"*^^ exceeds O(IOO), cohc !rcncc is lost and so oscillations are unobservable. Phase-averaging then 

1 

1 • 



sets (sin ) to zero and (sin^ MMjl^ iq i_ Implementing this averaging in Eq. (59) then leads to the classical 



expression 



(Pa.)pha«e-ave = ^(|A(n)|2) = ^ ^(An) 



11=1 



(n2 + v2)3 



sin^ 26»„ 



2 



(60) 



As with the quantum mechanical amplitudes in Eqs. (57) and (59), the bracketed expression pertains to the 
dNa weighting, while the analogous result for the model with dNp weighting is obtained by setting the brack- 
eted expression to {jS/ii^'^^Y . As we showed earlier, near a resonance the coherence length is increased by a factor of 
1/ sin 20. Thus, the phase-averaged Eq. (60) may not apply if LE lies near a resonance. 



SOME LIMITING CASES 



In this section we examine vacuum limits in 4D, and in our higher-dimensional 5D model. Vacuum limits are the 
expressions for oscillation probabilities that are valid at energies well below any resonance (s). The vacuum limits are 
different for the 5D model, and the standard 4D model. 
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A. The 4D Vacuum Limit 



The standard 4D vacuum expression for the two-state oscillation probability is the familiar one, 

2 



Pos(4D vacuum) = 



= sin^ 261 sin^ 



LAm'^ 
AE 



(61) 



The 5D Vacuum Limit 



At LE <^ {LE)^es (the notion of a resonant product {LE)^es will be discussed in section VI), smLSH/2 and sin^ 6 
given in Eqns. (9) and (11), respectively, assume standard 4D vacuum values. Thus, from Eqs. (59) and (57) we have 



Pas (5D vacuum) = Pas{4:D vacuum) x 



vn 



(„2 + ^2)3/2 



2 iSE 



(62) 



for the Gaufiian wo-distribution, and the bracketed expression replaced by (i/ii™^ for the Flat uo-distribution. 

Clearly, the 5D model is richer than the 4D model, and more complicated. Nevertheless, for consistency of our 
5D model, there must exist a limit in which the standard four- dimensional vacuum solution is obtained. Next we 
present this limit. 



C. The 4D Limit of the 5D Model 



Since in 4D all of spacetime comprises our brane, in the 4D limit the sterile neutrino, like the active neutrinos, is 
confined to the brane. Thus, the sterile neutrino exit angle ao = arctanu/n must be set to zero. This implies that the 
warp factor k must be set to zero. Equivalently, we may take |uo| to zero. According to Eq. (35), this again implies 
that k must go to zero. From Eq. (32) we see that e„ equals zero for all n in the A; = limit. The vanishing of e„ in 
turn implies that E'Res given by Eq. (58) is infinite, and so 5H and sin^ ^ given in Eqns. (9) and (11), respectively, 
assume standard 4D vacuum values. At this point, the 4D limit of the 5D probability has the form given in Eq. (62). 

Further simplification occurs. From Eq. (51) we sec that An is equal to unity in the fc = limit. Furthermore, the 
vanishing of e„ and Eq. (36) together reveal that e*'^'^''^") is an irrelevant, n-indcpcndcnt. miimodular constant. 

We are left to consider the mode integration / dN^, which must go to unity (one trajectory for the brane-bound 
sterile neutrino) in the limit A; — > 0. In the GauBian {D = G) case, we must take cr — > to reduce the p„-distribution 
to a delta function centered at p„ = 0, while in the Flat [D = F) case, we must take u^'^^ to reduce the 
Uo-distribution to a delta function at iio = 0. The zero limits of a or ii^'^^, and k, must be coordinated such that the 
mode integrals in Eqs. (44) and (46) go to unity in these limits. Let us define the ratios re = k/a and rp = k/ii^^. 
Prom Eqs. (44) and (46), we therefore require that 



1 I yf(^^^.-'(^ 



' , Gaufiian distribution 

(63) 

^ " I , Flat distribution. 



2 



Thus, a goes to zero as fc/r^, and u™*^ goes to zero as k/r*p, where r* is the solution of Eq. (63). For the Gaufiian 
case, Eq. (63) is transcendental. However, for the Flat case, the solution is simply r*p = 2/C,{2)(3L, where ((2) = 
X^n=i 1/?^^ — 7r^/6. The recipe for reduction of our model to the usual 4D vacuum solution is now complete. One 
takes fc ^ and r — > r*. 

Notice that n is an irrelevant parameter for motion on the brane. Accordingly, there is another, simpler, approach 
connecting our model to the 4D vacuum result. One may omit the sum on n and just set n equal to one. Then, the 
r* 's are solutions to the simpler normalization equations 



l(^§ELra'^e , Gaufiian distribution 

, Flat distribution . 



(64) 
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In particular, rp{n = 1) = 2/ [3L with this recipe. 

It should be noted that although our model in extra-dimensional warped space has a 4D vacuum limit, the extra- 
dimensional model is much richer in phenomenology. It includes two new physically-motivated parameters, the warp 
parameter k, and the wo-distribution parameter a or ii™'^^. One consequence of the dramatic difference is that there 
is no reason to expect the 4D vacuum parameters 6 and Am^ to correlate in value with the same parameters in our 
model. 



VI. THE NEAR ZONE 



Of particular interest is the "near zone'", which we define as 

~ ^ ) 01' equivalently, kL < 1 ["Near Zone"] . (65) 

The near zone is interesting for two related reasons. The first is that it is here in the near zone where an experiment 
can hope to see resonant behavior. From Eq. (58) we have that E^^^{n)/Am'^ = cos26'/2e„, so a resonant energy well 
above the sterile neutrino mass requires a small e„, as noted earlier in Eq. (8) and [5]. From the form of e{v/n) given 
in Eq. (34) or from Fig. 2 we have seen that small e„ means small v/n, exactly the condition employed to define the 
near zone. 

The second reason the near zone is interesting is that our formulas become especially simple in the near zone. The 
value of e„ at small v/n, given in Eq. (34) is v'^/Gn'^, which fixes the resonant energies, for E^^g{n) AmP, to be 



EResin) = — \/l2Am2 cos 26*. (66) 

k-L/ 

Something very interesting has emerged. The two parameters under experimental control are the energy E and the 
baseline L. We define the "resonant product of E and L" as 

1 



(L£)Res = T Vl2Am2 cos 26* . (67) 



Then we have that 

ElJn) nnLE)l^. ^""^ 

in the near zone. Substituting this result into the generalizations of Eqs. (11) and (9) to the geodesic modes, we 
obtain 



2E 

and 



A 777 2 / 

5H^ = —— 4 / sin^ 261 + cos^ 20 



{LEY 



n^{LE)l.. 



(69) 



. sin 26* 

sin26i„ = , . (70) 



sin^ 261 + cos2 261 



1 _ -ikMf 



It is now clear that in the near zone, where -E^es ^ Am^ or equivalently A;L <C 1, the location of the resonance 
depends not just on E but also on L, in the very simple combination LE. 
The amplitude for the n*'* mode contains the product of sines 

sm2&„ sm — - — — > sm ^^TrA/ -^sm26'j , (71) 

where the rhs is the maximum value of the product, attained whenever a resonant value of LE is reached, in which 

case the square root in Eq. (69) or (70) is minimized. Note that the resonant value is independent of n, and is bounded 
in magnitude by unity. Thus, this product cannot contribute significantly to the convergence of the sum, which we 
will discuss shortly. 

Further simplification of formulas appear in the near zone. The expansion of e„, Sn and An in powers of small 
v/n are already given in Eqs. (34), (37), and (51), respectively. Eqs. (47) and (51) yield the small v limits for the 
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Mo-distributions J dNo and J dNp. The ultimate result is the near zone limit for the amplitude of Eq. (57), which to 
leading nonzero order in v'^ is, for the case of the Gaufiian ttg-distribution. 



C30 



2 PE 



sin2^„ sin 



(72) 



where sin2^„ and sm{L5Hn/2) are given by the their near zone expressions in Eqs. (70) and (69), which show the 
explicit dependence on the resonant product (Li?)Res defined in Eq. (67). For the Flat UQ-distribution, the result is 
even simpler, as the bracketed expression in Eq. (72) is replaced by (i/u^^^ for this case. The near zone expression 
for An is 



An 



1, 



3/3n3 



(73) 



As discussed below Eq. (51), if coherent oscillations are sought, then v < 1/20 implies that min{- • • } in Eq. (73) is 
unity, and so An = 1 for all n. 

It is worth noting that the oscillation probability of Eq. (72) (and the mixing probability that results from this 



equation by setting (sin ^'^^" ) = 0, 
can be written as 



sin^ MMjl"^ = ^)) have an n-independent factor of v'^. Using Eq. (66), this factor 



3 Ato^ cos 29 



(74) 



The meaning 

.2 /p2 



of this exercise is that the oscillation probability (and mixing probability) is proportional to 
Am^/E^^^{n = 1). Thus, at least in the case of the Flat wo-distribution for which there are no other factors ex- 
ceeding unity, the energy of the first resonance cannot be too much larger than the sterile neutrino mass if the mixing 
is to be observable. For example, for an oscillation probability of 0.3%, the energy of the first resonance must occur 
at about thirty times the sterile neutrino mass. 



A. Convergence of the sum of amplitudes 

It is well known that the Riemann sum C(s) = X]^i(l/"'^) convergent for s > 1, and divergent for s < 1. In 
this subsection, we investigate the convergence properties of Pas as given in Eq. (72). For simplicity, we will assume 
the more common outcome of Eq. (73) and set An = 1. Since An is bounded from above by unity, this assumption 
makes little difference in what follows. Also, the two sine factors in Eq. (72) are bounded in magnitude by unity, and 
so will not be included in what follows. We are left to analyze the rates of convergence of 

N f a/^^ e~^(^) , Gaufiian uq- distribution 

^tSit , Flat Mo -distribution , 

for our two Mo-distributions. The analysis is quite different for the Flat and Gaufiian distributions. 



1. Convergence of Flat uo- distribution 

With the Flat Mo-distribution, it is clear that the summation is convergent, with each term in Eq. (75) weighted by 
1/n^. The convergence of the sum is modulated somewhat by the two sine factors in Eq. (72). 

It is true that any experiment will have a resonant mode, at nR ^ {LE) / {LE)jies, where {LE)^^^^ is the constant 
defined in Eq. (67). However, for the modes with n ^ n/f, the mixing factor sin20„ is very suppressed, while for the 
modes with n ^ ur, both sine factors approach the 4D vacuum values. Since the resonant amplitude is suppressed 
by i/nj^, the resonance feature will not be observable unless njj is a small integer. Thus, for the Flat UQ-distribution 
case, the first few modes will dominate the sum, and we may expect that the n = 1 resonance at LE = (LiJ)Res is 
the most manifest. 
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Experiment 


Baseline L [m] 


Energy E [MeV] 


Product LE [kmGoV] 


Xad [m] 


N[^°^ = L/X4D 


Bugey 


25 


1-6 


(2.5 - 15) X lO"'' 


2.5-15 


1.7-10 


KARMEN 


17.7 


20-52.8 


(0.35-0.93) xlO-^' 


50-130 


0.13-0.35 


LSND 


30 


20-50 


(0.6 - 1.5) X 10-3 


50-125 


0.24-0.6 


SNSosc (ORNL) - near 


18 


15-53 


(0.27 - 0.95) X 10-=^ 


38-130 


0.14-0.47 


SNSosc - far 


60 


15-53 


(0.90 - 3.2) X 10"=* 


38-133 


0.47-1.6 


FINcSSE (FNAL) near 


20 


500-1000 


(1 - 2) X 10"^ 


1.300-2500 


0.008-0.015 


FLN(-S.SE far 


100 


•lOO-lOOO 


().()")-(). 10 


i:500-2:)00 


0.0 10-0. 077 


MiniBooNE 


540 


300-800 


0.16-0.43 


740-2000 


0.3-0.7 


CDHS 


755 


> 1000 


> 0.76 


> 2500 


< 0.3 


Kamland (dominant reactor) 


180 X 10^ 


1-7 


0.18-1.26 


2.5-17.6 


(1 - 7) X 10* 


SuperK (subGeV atmos.) 


(10 - 10*) X 10=' 


500-1000 


5.0-10* 


1200-2500 


4-8000 


Minos 


735 X 10=^ 


(1 - 50) X 10^ 


730 - 3.7 X 10* 


(2.5 - 125) X 10* 


6-300 



TABLE I: Baseline L, energy range E, product LE, 4D vacuum oscillation length A413, and the number of vacuum 

oscillation lengths in the baseline, for several neutrino data sets [3, 4, 15-22], in order of increasing LE. For A4C and A'^^*^', a 
neutrino mass-squared difference of 1 eV^ is Eissumed. The SNSosc experiment is a neutrino detector proposed for the Spallation 
Neutron Source at Oak Ridge National Laboratory. 



2. Convergence of Gaufiian iio -distribution 



For the GauBian -uo-distribution, we may write 



n=l ^ ' 



_f l3EkL \ 



(76) 



Here the rate of convergence of the series depends critically on the slope of the GauBian e"^^/^. The magnitude of 
the slope monotonically increases from zero at 2; = and 00, to its maximum value of 1/y/e at 2; = 1 (the inflection 
point "IP"). As n increases, z decreases in inverse proportion. Thus, in the region z > 1, increasing n will drive the 
magnitude of the slope to larger values, and the (n + 1)** term in the series will be at least as large as 



(«+i) 



times 



the previous term. On the other hand, for 2; < 1, increasing n will drive the magnitude of the slope to smaller values, 
and the series will converge faster than 1/n. The consequence is that the series is mathematically convergent, but all 
amplitudes with > 1 will contribute significantly to the sum. An equivalent statement is that all amplitudes up to 



n ~ n/p = 



(3LEk 
2a 



1000/3 



LE 



0.4 kmCeV 



k 



10- 



Au\ 
urn) 



(77) 



contribute significantly to the sum. One sees that, depending on the variables LE, k and l/cr ~ Am, the number of 
significant amplitudes, n/p, may be large. Moreover, the number of significant amplitudes scales with LE. 



B. Near zone phenomenology 

It is illuminating to present examples of the new phenomenon arising from the L_E-resonancc of the warped extra- 
dimensional model. What becomes important for each experiment is the product of baseline times neutrino energy. In 
Table (I) we list some baseline and energy parameters for several recent and proposed experiments [27] . One possible 
phenomenology would be a resonant value of LE between the LSND and MiniBooNE values of 2.5 x 10"^ km GeV 
and 2.5 x 10"^ km GeV, respectively. With such a resonant value of LE, active-sterile vacuum oscillation, or even the 
resonance, could explain the LSND excess, while no observable active-sterile mixing would be expected in MiniBooNE, 
or in any other longer-Li? experiments. In particular, active-sterile mixing would be suppressed in the SuperK data 
sample of atmospheric neutrinos. 

We leave for future work the matching of the model with existing neutrino oscillation data. Here we explore the 
LE resonant phenomenon without reference to data. As with the discussion above of convergence rates for amplitude 
sums, the discussion of phenomenology bifurcates for the two Uo-distributions which we employ. 
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J as 



L[m] 

FIG. 3: Oscillation probability as a function of the baseline L for an assumed flat distribution of uq. Parameter values are 
^max = ^ ^ 1^ sin^26l = 0.003, k = 5/km, E = MeV, and Am^ = 2500 keV^ The resuhing value of (LJ5)Rea is 35 m MeV. 
For our choice of E, the resonance peaks occur at L = n {LE)ncB/ E — 35n m, n — 1, 2, 3 ... , with the principal resonance 
corresponding to n = 1. 



1. Phenomenology with Flat iio- distribution 

With the Flat ?io-distribution, the 1 /n^ falloff in amphtude imphes that the small values of n provide the dominant 
contributions to the oscillation probability. In particular, n — 1 defines the most dominant "principal resonance", as 
can be seen in Fig. 3. The peaks in the oscillation probability are observed at distances equal to integer multiples 
of {LE)yics/ E. The heights of the successive peaks decrease rapidly with increasing n, roughly as 1/n^. We have 
chosen parameters so that the peak of the principal resonance is of order of the mixing parameter inferred from the 
LSND excess. Attainment of this order of magnitude within the context of the new model is a demonstration of the 
potential relevance of the model to the real world. 

To the extent that the principal resonance is dominant, the phenomenology is relatively straightforward. At fixed 
baseline, the principal resonant mixing occurs at a unique energy, with 5D vacuum oscillations pertaining at energies 
well below this principal resonant energy, and suppression of mixing pertaining at energies well above. This part 
is standard oscillation phenomenology [28]. What is truly new is that at fixed neutrino energy, resonant mixing 
occurs at one baseline value L, with vacuum oscillations pertaining at much smaller baselines and suppression of 
mixing pertaining for much larger baselines. So the sterile neutrino decouples from active neutrino physics at either 
sufficiently large E or sufficiently large L. The principal resonance is not an energy resonance or a length resonance, 
but rather a product LE resonance. So a more careful statement is that once the resonant value of LE has been 
identified experimentally, then 5D vacuum oscillations apply at smaller LE, and active-sterile mixing is suppressed 
at larger LE. 

One consequence of the principal resonance condition LE « (LE)^^cs is that, should an experiment observe this 
resonance, then in the context of this model one may infer the value of the warp-parameter. From Eq. (67), we have 
that k = Vl2Am^cos20/ {LE)^s- Corroborative information may be carried by the height of the principal resonance 
peak. From Eq. (67), one easily derives that on the principal resonance, {vE)^^s — VSAm^ cos26', and so 

(-Pas)first peak = [An{n = 1)]^ 

The L and E appearing in this formula are correlated by the resonance condition Eq. (67). 



kL 



3Am2 008 26* 



sm 



/LAm^ sin 26* 



AE 



(78) 



2. Phenomenology with Gaufiian uq- distribution 

The oscillation probability versus baseline for a GauBian UQ-distribution is shown in Fig. 4 (red curve). With the 
Gaufiian UQ-distribution, there may be many modes making significant contributions to the overall amplitude, in 
which case an analytical analysis seems difficult. The number of significant modes is indicated in Eq. (77). For the 
parameters we use in Fig. 4, the number of significant modes is of the order of njp = 27 (L/200m). For the principal 
resonance at L = 37 m, n/p ~ 5, whereas for the fifth resonance at L = 5 x 37 m, one gets njp 25. 
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FIG. 4: Oscillation probability as a function of the experimental baseline, for the Gaufiian distribution (red and green curves). 
The green curve presents the phase-averaged oscillation probability, and the sinusoidal blue curve presents the probability 
as given by the standard 4D vacuum formula for oscillations between sterile and active neutrinos. Parameter choices are 
sin^26l = 0.003, k = 5/(10** m), £ = 15 MeV, Am^ = 64 eV^ and a = 100 eV. The resulting value of {LE)nes is 550 m MeV. 
For our choice of E, the resonance peaks are found at the multiples L = n{LE)-R.es/ E = 37n m, n = 1, 2, 3 • • • , with the principal 
resonance corresponding to n = 1. 

Still, Fig. 4 reveals an oscillation probability with peak heights and peak structures not dissimilar to the case of the 
flat distribution (although parameter choices are quite different in the two cases). The peaks are again found at 
multiple integers of the principal resonant product {LE)iies/ E. 

Shown in green in Fig. 4 is the phase-averaged oscillation probability. For our parameter choices, the phase-averaged 
plot provides a faithful representation of the oscillation probability even for the low-energy principal resonance. 

We also show (in blue) in Fig. 4 the 4D vacuum model of active-sterile oscillation, given in Eq. (61). It is well-known 
that the two-state active-sterile 4D vacuum solution, while it can explain the LSND data in isolation, cannot explain 
the comprehensive sets of data, including null results of other short-baseline experiments [14]. Still, we show it for 
comparison to our extra-dimensional model. Large differences between the two models are apparent. In particular, 
the 4D vacuum result has maxima where the argument of sin^(LAm^/4£^) equals odd multiples of 7r/2, i.e., at 
L = {2ttE / Am'^) x q, with q = 1,3,5,.... On the other hand, the extra-dimensional model has maxima roughly 
where sin0 equals unity, i.e., at LE = {LE)b^cs x n, with n — 1, 2, 3, ... . Hence, within a fixed baseline or energy 
range, there are on average twice as many peaks in the extra-dimensional model as in the 4D vacuum model. 

The novelty of the LE resonance is best illustrated in a three-dimensional plot, with oscillation probabilities rising 
as topography on a two-dimensional L-E plane. Contours of constant LE appear on this plot as hyperbolas symmetric 
about the diagonal line L = E. With the Gaufiian WQ-distribution, inspection of the near-zone formula in Eq. (72) 
reveals that scaling in LE (recall that v = kL/2) is broken only by the explicit L in the argument of sm{LSHn/2) (or 
equivalently, by an explicit E if we write L = [LE)/E). Thus, we expect the oscillation probability to roughly peak 
at resonant values of LE, i.e., along hyperbolic contours of constant LE. Such peaking is clearly evident in Fig. 5. 
The hyperbolic contours of resonant probability display a marked contrast from, say, MSW resonances which occur 
at resonant energy values independent of baseline L; and the hyperbolic iso- LE contours contrast with 4D vacuum 
maxima, which occur at sin(LAm^/4i?) = 1, i.e., at rays in the L-E plane given by L/E = q2TT/ Am?, q = 1, 3, 5 . . . 

We end this section with a comment on the possible slow convergence of the mode sum in the case of the Gaufiian 
MQ-distribution. If the parameters are such that njp = (3LEk/2cr defined in Eq. (77) is large compared to one, then 
there are several or more contributing modes, and the concept of principal resonance gives way to a contributing band 
of resonances about the mode number 



with {LE)-iics given in Eq. (67). For the resonances in this band, the amplitudes are enhanced compared to their 4D 
vacuum analog by the factor 



LE 



(79) 



nRband = 



{LE)ncs 




(80) 
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FIG. 5: Oscillation probability (vertical) in the L-E plane, with the same parameters as in Fig. 4. Units of L and E are m and 
eV, respectively. 

For small 9, but large ttN^^"^ sin 26* > 1, this enhancement of the ampUtudes in the resonant band can be quite large, 
typically of order 1/26*. This enhancement is presumably operative in Figs. 4 and 5, where, as we have stated above, 
nip ranges from 5 at the first resonance band to 25 at the fifth resonance band. 

VII. DISCUSSION AND CONCLUSION 

The basic feature of the new model presented herein is that longer neutrino travel times associated with longer 
baselines on the brane allow the off-brane geodesic of the sterile neutrino to plunge deeper into the bulk and experience 
a greater warp factor. As a consequence, the shortcut parameter e = St/t for the sterile neutrino increases with the 
baseline, corresponding in turn to a decreasing resonance energy in the effective Hamiltonian of the two-neutrino 
system. In addition, there are higher energy/longer baseline resonances resulting from additional classical geodesies. 
We solved the geodesic equations and identified the countably infinite number geodesies giving rise to the countably 
infinite number of resonances. Then we performed a path-integral-weighted sum over the amplitudes resulting from 
each of the geodesies. 

In the "near zone", defined as baselines short on the scale of the warp factor fc~^, the resonance condition is that 
the product of baseline and energy, LE, be an integer multiple of a fundamental value (Li?)Res ~ rrig/k. That the 
brane-bulk resonances in the "near zone" expansion depend on the product of the energy and baseline, rather than 
on the energy alone as with the MSW matter-resonance, is a novel feature of our model. 

Whether or not the higher resonances coming from the additional geodesies contribute significantly depends on 
the initial distribution of sterile velocities wq transverse to the brane. We considered two very different distributions 
for initial -uq, a- GauBian distribution with a width related via the uncertainty principle to the brane width, and a 
Flat distribution as might arise in a path-integral approach. Both approaches led to qualitatively similar oscillation 
probabilities, although the parameters in the two approaches, and the physical interpretation of some of the param- 
eters, were quite different. Also, the contributions of higher resonances are more suppressed in general for the Flat 
distribution than for the Gaufiian distribution. 

Since higher-ii? resonances are suppressed, and active-sterile neutrino mixing is suppressed for LE above the 
resonant values, sterile neutrinos taking shortcuts in the extra-dimensional bulk decouple from active neutrinos in 
long-baseline experiments. Thus, no active-sterile mixing is expected in atmospheric data, in MINOS or CDHS. All 
explanations proposed so-far for the LSND and MiniBooNE anomalies assumed baseline-independent oscillations and 
mixing. Difficulties accommodating longer-baseline data were encountered in these models. These difficulties do not 
immediately extend to scenarios with warped extra dimensions, as developed here. In fact, the failure of previous 
models to reconcile short baseline data such as LSND with longer baseline data might be construed as favoring the 
extra-dimensional shortcut scenario. Finally, we mention that the bulk shortcut scenario might even relieve some of 
the remaining tension between the LSND and KARMEN experiments. Since LSND has almost twice the baseline of 
KARMEN, the bulk shortcut model opens more parameter space for accommodating the two experiments. 

In conclusion, we have found that in scenarios with sterile neutrinos taking shortcuts in the extra-dimensional bulk, 
the shortcut is baseline L dependent as well as energy E dependent. Resonances occur as a function of L as well E, 
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and may lead to neutrino anomalies. Above the contributing resonances, i.e., at large L or large E, the mixing sterile 
neutrinos with active neutrinos is suppressed. In the context of this model, existing data on neutrino oscillations and 
anomalies need to be reanalyzed. We intend to compare the model to the world's data in future work. 
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